State-insensitive trapping and guiding of cesium atoms using a two-color evanescent 

field around a subwavelength-diameter fiber 
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We calculate the optical potentials, i.e. the light shifts, of the ground and excited states of atomic 
cesium in a two-color evanescent field around a subwavelength-diameter fiber. We show that the 
light shifts of the 65*1/2 <-* 6P3/2 transitions can be minimized by tuning one trapping light to around 
934.5 nm in wavelength (central red-detuned magic wavelength) and the other light to around 685.5 
nm in wavelength (central blue-detuned magic wavelength). The simultaneous use of the red- and 
blue-detuned magic wavelengths allows state-insensitive two-color trapping and guiding of cesium 
atoms along the thin fiber. Our results can be used to efficiently load a two-color dipole trap by 
cesium atoms from a magneto-optical trap and to perform continuous observations. 

PACS numbers: 32.80.Pj,32.80.Lg,03.75.Be,03.65.Ge 



I. INTRODUCTION 

There is currently a strong interest for the manipula- 
tion of individual neutral atoms in microscopic (at sub- 
wavelength size) optical dipole traps |l|,|2j]. The ability to 
manipulate atoms individually may open a way to con- 
trolled engineering of the quantum states of small sets of 
trapped particles, in order to encode and process infor- 
mation at the quantum level A red-detuned optical 
dipole trap works as follows: The ac Stark shift induced 
by a red-detuned trapping light lowers the ground-state 
energy of the atom proportionally to the local light in- 
tensity 4s] . The spatial dependence of the atomic poten- 
tial energy is therefore equivalent to the spatial depen- 
dence of the light intensity. The atom has the lowest 
potential energy in the region of maximum intensity and 
can therefore be trapped there. For very large detun- 
ing, typically several nanometers, the photon scattering 
rate becomes so low that the optical trapping potential 
is truly conservative. The optical dipole trap is usually 
loaded by precooled atoms from a magneto-optical trap 
(MOT). There are two schemes of microscopic optical 
dipole traps to store individual atoms. The first one is 
based on a strongly focused single Gaussian laser beam 
0. The second one is based on the use of a standing 
wave of two counterpropagating laser beams with equal 
intensities and optical frequencies 0- In both schemes, 
the size of the trapping region is less than the light wave- 
length. 

Recently, we proposed a new method for microscopic 
trapping of individual atoms 0, [fj . The method is based 
on the use of a subwavelength-diameter silica fiber with a 
single (red-detuned) light beam |{| or with two (red- and 
blue-detuned) light beams launched into it. The light 
waves decay away from the fiber wall and produce opti- 
cal potentials for neutral atoms. The red-detuned light 
wave produces an attractive potential. To sustain stable 
trapping, the atoms must be kept away from the fiber 
wall. This can be achieved by the centrifugal potential 
barrier in the one-color scheme Q or by the repulsive op- 



tical potential from the blue-detuned laser beam in the 
two-color scheme . The atom trapping and guiding 
occur outside the fiber. We have shown that stable trap- 
ping and guiding can be achieved only when the fiber 
diameter is smaller than the light wavelength. The great 
advantages of our scheme are (a) localization of atoms 
to a subwavelength region, (b) high efficiency to detect 
individual atoms, (c) high accessibility to the trapped 
atoms, and (d) achievement of strong coupling between 
light and matter 

Due to the conservative character of optical potentials, 
the loading of atoms into dipole traps requires the use of 
friction forces, which can be provided by the Dopplcr 
cooling mechanism in MOTs 0, The loading of a 
dipole t rap from a MOT is a dynamical process rich in 
physics [nj. A major obstacle to the efficient loading of a 
dipole trap is that the excited electronic states generally 
experience positive ac Stark shifts of comparable mag- 
nitude to the negative shift of the ground state. These 
light shifts with opposite signs change the resonant fre- 
quencies of the atoms in the dipole trap. When the spa- 
tial gradient of the trapping fields is high, the resonant 
frequencies of the atoms vary substantially with their po- 
sitions within the dipole trap. This effectively introduces 
a strong spatial dependence of the detuning between the 
atomic cooling transition and the MOT fields. The de- 
tuning may become large or positive. This prevents the 
efficient Doppler cooling. As a result, the Doppler cool- 
ing is generally incompatible with the dipole trapping. 
The extremely small volume of a microscopic dipole trap 
makes the problem even worse: the time-averaged num- 
ber of atoms in a microscopic dipole trap (loaded from a 
MOT) can be less than one. 

Katori et al. have proposed and demonstrated a trap 
loading scheme that helps overcome the above obstacle 
The idea is to use an appropriate wavelength of 
the trapping laser light (called the magic wavelength) at 
which (due to the specific multilevel structure of a real 
atom) the Stark shifts of the lower and upper levels of 
the cooling transition have the same value and the same 
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sign. It has been demonstrated for strontium atoms that 
the choice of a magic wavelength for a red-detuned far- 
off-resonance trap allows simultaneous Doppler cooling 
and dipole trapping and hence enables high loading effi- 
ciency of magneto-optically trapped atoms into the opti- 
cal dipole trap [Tl|. 

The spatial dependence of the atomic resonant fre- 
quencies also leads to additional complications in the 
monitoring (probing) processes. Kimble et al. have re- 
cently demonstrated for cesium atoms that the choice of 
a magic wavelength for a red-detuned far-off-resonance 
trap allows state-insensitive trapping and continuous ob- 
servation of trapped atoms . 

In this paper, we address the problem of minimizing 
the spatial dependence of the light shifts of the atomic 
transitions in a two-color dipole trap. We show that the 
light shifts of the 6Si/ 2 6P3/2 cesium transitions can 
be minimized by tuning one trapping light to around 
934.5 nm in wavelength (central red-detuned magic wave- 
length) and the other light to around 685.5 nm in wave- 
length (central blue-detuned magic wavelength). We 
calculate the optical potentials of cesium atoms in the 
ground and excited states in a two-color evanescent field 
around a subwavelength-diameter fiber. We show the 
possibility of state-insensitive two-color trapping and 
guiding. 

Before we proceed, we note that, due to recent de- 
velopments in taper fiber technology, thin fibers can be 
produced with diameters down to 50 nm [l3l IT^j . Thin 
fiber structures can be used as building blocks in future 
atom and photonic micro- and nano-devices. 

The paper is organized as follows. In Sec. [H] we re- 
view the general theory and present the basic equa- 
tions for light shifts and polarizabilities of atomic states. 
In Sec. IIIII we calculate the dynamic polarizabilities of 
the 65*1/2 ground state and the 6P 3 / 2 excited state of 
atomic cesium. In Sec. IIVI we calculate the light shifts 
of these states in a two-color evanescent field around a 
subwavelength-diameter fiber. Our conclusions are given 
in Sec. 



II. GENERAL THEORY OF LIGHT SHIFTS OF 
ATOMIC HYPERFINE LEVELS 



We consider the interaction between an atom in a fine- 
structure state \n) = \nLj) and an external electric field 
E. The combined Hamiltonian of the hyperfine interac- 
tion and the Stark effect is 



H — Vhf s + Vee, 



(1) 



where the operator Vhf s describes the hyperfine structure 
and the operator Vee describes the field-induced shifts 
of energy levels. 

The hyperfine interaction operator Vhfs is given by |l5| 



V B A T T I B H 6(1 ' J)2 ± 31 • J 2/(7 ± 1)J( J ± 1} 
V Ms HAI-J+hB —————— . 

(2) 

Here J is the operator for the total electronic angular 
momentum, I is the operator for the nuclear spin, and 
A and B are the hyperfine-structure (hfs) constants. For 
the level 6P3/2 of cesium (I = 7/2 and J = 3/2), we have 
@ A/2n = 50.34 MHz and B/2tt = -0.38 MHz. 

The Stark operator Vee is, in the second-order pertur- 
bation theory, given by [la ] 



Vee 



~a E 2 - ^a 2 QE 2 . 



(3) 



Here ckq and a 2 are the scalar and tensor polarizabilities, 
respectively. The scalar polarizability «o shifts all hyper- 
fine and magnetic sublevels equally. The tensor polariz- 
ability a 2 mixes the hyperfine and magnetic sublevels 
through the operator [l(| 



Q = 



3(u ■ J) 2 — J(J + 1) 
J(2J- 1) 



(4) 



Here u = ~E/E is the unit vector in the field direction. 

Due to the hfs interaction, the total electronic angular 
momentum J is not conserved. However, in the absence 
of the field, the total angular momentum of the atom, 
described by the operator F = J + 1, is conserved. In the 
basis of hfs states \FMp), the operator Vhf s is diagonal. 
Its nonzero matrix elements are 



(FMp I Vhfs I FMp ) 



1 



-KAK 



%K(K + 1)-2I(I+1)J(J+1) 
21(21 - 1)2 J(2J - 1) 



(5) 



where K = F(F + 1) - I {I + 1) - J(J + 1). 

The matrix elements of Q between two hfs states are 
given by (l6j| 



(FM F \Q\F'M' F ) 



(J+l)(2J + l)(2J + 3) 



J(2J- 1) 



1/2 



1 2 1 
fi -q n' 



x (-i) /+J+f - F '- MF VWTT)(2F 7 TT) 
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iUy)/y/2, 


U = 


u z 


and 



(6) 



Here u_i = (u x 

iu y ) I a/2 are the spherical tensor components of the field- 
direction vector u. The three-j symbols in eq. © require 
the two conditions q = /x + // = M' F — Mp. The second 
three-j symbol and the six-j symbol in eq. JBJ make Q 
not diagonal in F. 

Equations Q and © are valid for an arbitrary orien- 
tation of the electric field. In a particular case where the 
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electric field is aligned along the quantization axis z, i.e., 
E = Ez, eq. reduces to the form 



Q 



3 J 2 - J( J + 1) 
J(2J - 1) : 



(7) 



which is diagonal in J and Mj. Most of the previous 
works on the Stark effect in atomic excited states were 

devoted to this case ei nam. 

In the absence of the hfs interaction, the Stark shift of 
a fine-structure magnetic sublevel | JMj) can be written 
as AEm,j — ~(l/2)a(Mj)E 2 , with the polarizability 



a(Mj) = a + a 2 



3Mj - J (J + 1) 
J(2 J — I) ' 



(8) 



±3/2) 



In particular, for J = 3/2, we have a{Mj 
ao + a 2 and a(Mj — ±1/2) — ao — a 2 . 

The presence of the hfs interaction dictates the use 
of the hfs basis {\FMf)}- In this basis, the particular 
special form (JJJ of the operator Q is diagonal in Mp but 
not in F. The nonzero matrix elements of this operator 
are 



(FM F \Q\F'M F ) 



\I+J+F-F'-M F 



(J + I)(2J+1)(2J±3) 
J(2J- 1) 

V(2F+1)(2F' + 1) 



F 


2 


F> \ 




2 




M F 





—Mp ) 




I 





1/2 



(9) 



For each fixed value of Mf, there is a matrix with rows 
and columns labeled by F and F' , respectively. In par- 
ticular, for Mf = ±-F max , where F max = J + J, the ma- 
trix (FMf|Q|F'Mf) reduces to f x f. Hence, for the 
hfs states with the maximum values of F and |Mf|, i.e. 
the states \F = F max ,Mp = ±F max ), the Stark shift is 
A£?M max = -(l/2)aM m „-B 2 , where aM m „ = «o + ol 2 . 

In this paper we study the case where the field in- 
teracting with the atom is an optical field, that is, 
E = (Se- iut + £*e iut )/2. Here £ is the complex en- 
velope vector of the electric component of the light field. 
In this case, we have to remove fast optical oscillations 
from the right-hand side of eq. J2J) by averaging it over 
an optical period. This procedure leads to the following 
expression for the operator of the dynamic Stark effect: 



V EE = - \a \£\ 2 - \a 2 £ \FM F )(F'M' F \ 



qflfl' 

T5 



FM F F'M' F 

1 2 I 

i" -q p! 



(J + l)(2J+l)(2J + 3) 
J(2J- 1) 



1/2 



(10) 



X ^y+J+F-F'-Mr ^( 2F + j^p + ^ 
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Here £_i = (£ x - i£ y )/y/2, £ Q = £ z , and £i = -(£ x + 
i£ y )/^/2 are the spherical tensor components of the field 
envelope vector £. In general, the Stark operator Vee is 
not diagonal in the hfs basis {\FM F }}, and consequently 
neither is the Hamiltonian . To find the shifts of the 
hfs sublevels, we must diagonalize this Hamiltonian. The 
Stark shift induced by an optical field is called the ac 
Stark shift or the light shift. 

In the case of optical fields, in addition to the time av- 
eraging procedure, we also have to use the dynamic polar- 
izability instead of the static one. The static scalar and 
tensor polarizabilities of an atomic fine-structure state 
have been derived systematically |17|. The microscopic 
expression for the linear susceptibility, which is related to 
the dynamic polarizability, is well known in the literature 
[T^ | . Combining the previous results , we write the 

dynamic scalar polarizability ao and the dynamic tensor 
polarizability a 2 of an atomic fine-structure state |n) as 



2/ft \- , ||n|| 
3(27+1) V " " 



v ^n'»K'»"^ 2 +7n'n/4) , , 

(uj 2 , -UJ 2 +1 2 , /4) 2 ± 7 2 , lu 2 1 ' 



and 



a 2 



5J(2J- 1) 



h V6(J+l)(2J±l)(2J + 3) 
[J+ J> \ J 1 J' 



1/2 



x 2J("1) 



1 J 2 

^n'«K 2 t ' ?t -^ 2 ±7n'n/ 4 ) 

(u 2 , -w 2 ±7 2 , /4) 2 ±7 2 , uj 

\ n'n 1 in'nl t 1 in'n 



n'\\D\\n) 2 



(12) 



respectively. Here uj is the frequency of the light field, 
(n'||D||n) is the reduced electric dipole matrix element 
for the transition between the fine-structure states \n') 
and \n), oj n ' n — u> n i -w„ is the transition frequency, and 
In'n — 7n' + In is the linewidth (twice the dephasing 
rate) and is given as the sum of the population decay 
rates j n i and j n of the states. We note that, when uj is 
several linewidths off resonance with the corresponding 
transition, the effect of 7„' n in the above formulae can 
be neglected. We also note that, for the ground states 
of alkali-metal atoms, which correspond to J = 1/2, the 
tensor polarizability is vanishing, that is, a 2 = 0. 



III. NUMERICAL RESULTS FOR CESIUM 
ATOMS IN A LINEARLY POLARIZED 
PLANE- WAVE LIGHT FIELD 

A. Dynamic polarizabilities of the 6S1/2 ground 
state and the 6P3/2 excited state 

The polarizabilities of the ground and excited states 
of atomic cesium have been calculated in a large number 
of works. However, most of the previous calculations 
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were devoted to the static limit EEE3. Recently, in 
order to search for a red-detuned magic wavelength for a 
far-off-resonance trap, the light shifts of the ground and 
excited states of atomic cesium have been calculated ■ 
However, the results for the dynamic polarizabilities have 
not been explicitly provided. 

To search for red- and blue-detuned magic wave- 
lengths, we first calculate the polarizabilities of the 
ground and excited states of cesium as functions of the 
light wavelength A, using eqs. i|II|) and (|12|l . The 
calculations for the polarizability of the 6Si/ 2 ground 
state incorporate the couplings 6S1/2 «-» (6-ll)Pi/2,3/2- 
The calculations for the scalar and tensor polarizabili- 
ties of the 6P3/2 excited state incorporate the couplings 
6P3/2 <-> (6-15)Si/ 2 and 6P 3 / 2 <-> (5~ll)As/ 2 ,5/2- Rel- 
evant parameters are taken from a number of sources 

HI IHI2II23. 

We plot in Fig. ^the dynamic polarizability a — ao of 
the ground state 6Si/ 2 . As seen, in the region A > 500 
nm, the profile of a has two closely positioned resonances, 
corresponding to the transitions between 6 Si/ 2 and 6P1/2 
(Di line, wavelength 894 nm) and between 6Si/ 2 and 
6-P3/2 (Pi line, wavelength 852 nm). The effects of the 
other transitions are not substantial in this wavelength 
region. 
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FIG. 1: Polarizability a of the ground state 6S1/2 in atomic 
cesium as a function of the light wavelength A. 

We plot in Fig. [2 the scalar polarizability ao and the 
tensor polarizability ai for the excited state 6P 3 / 2 . The 
figure shows that both ao and a2 have multiple reso- 
nances in the region A > 500 nm. The most dominant 
resonances are due to the transitions from 6P3/2 to (6- 
8)S 1/2 and (5-8)^3/2,5/2- 



B. Blue- and red-detuned magic wavelengths 

We now search for magic wavelengths, at which the po- 
larizabilities and consequently the light shifts of the rele- 
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FIG. 2: Scalar polarizability ao (a) and tensor polarizability 
«2 (b) of the excited state 6P 3 / 2 in atomic cesium as functions 
of the light wavelength A. 



vant upper and lower states are almost equal, leading to 
the minimization of the shift of the atomic transition fre- 
quency For this purpose, we plot in Fig. [3]the sum 
(a) and difference a^—a-i (b) of the scalar and ten- 
sor polarizabilities of the 6P3/2 excited state (solid lines) 
together with the polarizability of the 6Si/ 2 ground state 
(dashed lines). As mentioned in the previous section, the 
quantities ao + a2 and ao — a2 are the polarizabilities of 
the fine-structure magnetic sublevels with Mj — ±3/2 
and Mj — ±1/2, respectively, in the case where the field 
is linearly polarized along the z axis. Therefore, although 
the total polarizability of the 6P3/2 excited state is a ten- 
sor, the quantities ao + a2 and ao — a2 characterize the 
boundary magnitudes of the total polarizability. 

As seen from Fig. [21 the sum ao + a2 and difference 
ao — a2 for the 6P3/2 state cross the polarizability of the 
6S1/2 state at slightly differing wavelengths of 938 nm 
and 931 nm, respectively. These crossing points are red- 
detuned from the D\ and D 2 resonance lines. They are 
spread around the central red-detuned magic wavelength 
Xr = 934.5 nm, in agreement with the recent result of 
McKeever et al. for atomic cesium [T^ . 

We observe in Fig. EJthat, in addition to the crossings 
of the ground- and excited-state polarizabilities on the 
red side of detuning, there are several crossings on the 
blue side. The two blue-detuned crossings that are clos- 
est to the Di and D2 resonance lines occur, in the case 
of ao±a2, at the wavelengths of 684 nm and 612 nm [see 
Fig. Ola)] and, in the case of oq — at the wavelengths 
of 687 nm and 614 nm [see Fig. E£b)]. The difference be- 
tween the positions of the ao + a2 and ao — a2 crossings is 
rather small. The first blue-detuned crossings are spread 
around the central magic wavelength As = 685.5 nm. 
The second blue-detuned crossings are spread around the 
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FIG. 3: Comparison between the polarizabilities of the excited and ground states of atomic cesium. The sum (a) and difference 
(b) of the scalar and tensor polarizabilities of the 6P3/2 excited state are shown as functions of the light wavelength by the 
solid lines. The polarizability of the 6S1/2 ground state is shown by the dashed lines. 



central magic wavelength X' B = 613 nm. Thus we can 
minimize the light shifts of the 65*1/2 6P3/2 transi- 
tions of cesium atoms in a blue-detuned light field by 
tuning the light field to around an average wavelength 
Xb = 685.5 nm or X' B — 613 nm. Concerning the problem 
of trapping and guiding atoms around a subwavelength- 
diameter fiber, the first blue-detuned magic wavelength 
Xb = 685.5 nm is more favorable than the second blue- 
detuned magic wavelength X' B = 613 nm. One of the 
reasons is that the first wavelength is closer to the D\ 
and D2 resonance lines and hence leads to a larger cou- 
pling strength. In addition, the first wavelength satisfies 
better the single-mode fiber condition and gives a longer 
evanescent-wave penetration length. Therefore, we focus 
on the first blue-detuned magic wavelength but not on 
the second one. 

We show in Figs. Ufa) and^Jb) the polarizabilities of 
the 65i/ 2 and 6P 3 / 2 states in the vicinities of the cen- 
tral blue-detuned magic wavelength Xb = 685.5 nm and 
the central red-detuned magic wavelength A^ = 934.5 
nm, respectively. As seen, around Xb and A^, the po- 
larizabilities of the ground and excited states cross each 



other. The signs of the polarizabilities in the vicinities of 
Xb and Xr are negative and positive, respectively. The 
magnitudes of the polarizabilities in the vicinities of Xb 
and Xr are on the order of —600 a.u. and 3000 a.u., re- 
spectively. The magnitudes of the polarizabilities in the 
vicinity of Xb are about five times smaller than in the 
vicinity of Xr. The reason is that As is farther from the 
Di and D 2 resonance lines than Xr. 



C. Light shifts due to a linearly polarized light 

We calculate the light shifts of the transitions from the 
6P 3/2 FM F sublevels to the 6S 1/2 F'M' F sublevels. The 
light shift of an atomic transition is the difference be- 
tween the light shifts of the upper and lower levels. The 
light shifts of the hfs sublevels \FMp) of the 6P 3 / 2 excited 
state are determined by diagonalizing the Hamiltonian 
(JIJ , which includes the hfs energy and the Stark inter- 
action energy (|10|l . The Stark energy of the 6P 3 / 2 state 
is produced by the scalar polarizability cvo(6-p3/2) and 
the tensor polarizability 02(6/3/2)- All the hfs sublevels 
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FIG. 4: Blue-detuned (a) and red-detuned (b) magic wave- 
lengths for the 6S1/2 <-» 6P3/2 cesium transition. The sum 
and difference of the scalar and tensor polarizabilities of the 
6P3/2 excited state are shown by the thick and thin solid lines, 
respectively. The polarizability of the 651/2 ground state is 
shown by the dashed lines. 



FIG. 6: Light shifts of the transitions from the 6P3/2FMF 
sublevels to the 6S1/2F' M' F sublevels in cesium vs. the field 
intensity. The values of the wavelength of the field are chosen 
to be A = 684 nm (a) and A = 938 nm (b). The electric com- 
ponent of the field is linearly polarized along the quantization 
axis z. We show only the results for the transitions involving 
the excited-state sublevels that are split from F — 5. 




X [nm] 



FIG. 5: Light shifts of the transitions from the 6P 3 / 2 FMf 
sublevels to the 6S1/2F' M' F sublevels in cesium as functions 
of the light wavelength in the vicinities of the central blue- 
detuned magic wavelength (a) and the central red-detuned 
magic wavelength (b). The electric component of the field is 
linearly polarized along the quantization axis 2. The intensity 
of the field is 1 MW/cm 2 . We show only the results for the 
transitions involving the excited-state sublevels that are split 
from F = 5. 



\F'M' F ) of the 6S l 1 / 2 ground state have the same light 
shift, produced by the scalar polarizability ao(65i/2). 
For simplicity, we assume in this subsection that the elec- 
tric component of the light field is linearly polarized along 
the quantization axis z. In addition, we limit ourselves 
to the transitions from the excited-state sublevels that 
are split from F — 5. 

In Fig. [S] we plot the shifts of the transition frequen- 



cies as functions of the light wavelength in the vicinities of 
the central blue-detuned magic wavelength Xb = 685.5 
nm (a) and the central red-detuned magic wavelength 
A_r = 934.5 nm (b). As seen, the light shifts cross zero 
at around Xb and A^, with positive slopes. At Xb and 
Xr, the light shifts range from —4 MHz to 3.1 MHz and 
from -21.2 MHz to 26.3 MHz, respectively. The range 
of the light shifts in the vicinity of Xb is several times 
smaller than that in the vicinity of Xr. This is due to 
the difference between the magnitudes of the polarizabil- 
ities around Xb and Xr. 

In Fig. H3 we plot the shifts of the transition frequen- 
cies as functions of the field intensity. We choose the 
wavelengths A = 684 nm (a) and A = 938 nm (b) for the 
field. These values are close to the exact values of the 
blue- and red-detuned magic wavelengths for the tran- 
sitions 6P 3/2 F = 5M F = ±5 <-> 6S 1/2 F' M' F (with the 
maximum values of F and |Mp|). Our numerical cal- 
culations show that the light shifts of these transitions 
are indeed small. They are less than 2 MHz even when 
the field intensity is as high as 3 MW/cm 2 . A more pre- 
cise tuning can, in principle, reduce the light shifts of 
these transitions to zero. The light shifts of the other 
transitions are more substantial but not very large. At 
an intensity of 3 MW/cm 2 , the light shifts of the tran- 
sitions from the F = 5 Mf = sublevel are —23 MHz 
in the case of Fig. [Bfa) and 146 MHz in the case of Fig. 
HJb). In general, it is possible to individually minimize 
the light shifts of the transitions from the upper sublevels 
6P3/2FMF with arbitrary fixed values of F and Mp us- 
ing appropriate choices of magic wavelengths. 
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IV. LIGHT SHIFTS OF THE TRANSITIONS OF 
ATOMIC CESIUM IN A TWO-COLOR 
EVANESCENT FIELD AROUND A 
SUBWAVELENGTH-DIAMETER FIBER 

Consider a cesium atom moving outside a thin single- 
mode optical fiber that has a cylindrical silica core of 
radius a and refractive index n\ and an infinite vacuum 
clad of refractive index n% = 1. To produce an optical po- 
tential with a trapping minimum sufficiently far from the 
fiber surface, we use two laser beams propagating along 
the fiber in the fundamental modes 1 and 2 with differing 
frequencies u>i and u> 2 , respectively (with wavelengths Ai 
and A2 , respectively, and free-space wave numbers fci and 
&2j respectively) ||. To make the potential cylindrically 
symmetric, the laser beams are circularly polarized at 
the input. In the vicinity of the fiber surface, the polar- 
ization of the transverse component of each propagating 
field rotates elliptically in time, the orbit rotates circu- 
larly in space, and the spatial distribution of the field 
intensity is cylindrically symmetric [? ]. For certainty, 
we assume that circulation of photons around the fiber 
axis z is clockwise. 

Outside the fiber, in the cylindrical coordinates 
{r, ip, z}, the cylindrical components of the envelope vec- 
tor £ of the electric field in a fundamental mode with 
clockwisely rotating (circulating) polarization are given 

by [? ] 

£ r = tA[{l- S )K (qr) + {l + s)K 2 (qr)}e^ z - ip \ 
£ v = A[(\- S )K a (qr)-(\ + s)K 2 (qr)]e^ z ^\ 

£ z = A^K^qrY^-^ . (13) 
P 

Here (3 is the longitudinal propagation constant deter- 
mined by the eigenvalue equation for the fiber mode 
with the free-space wavenumber k = lu/c, the param- 
eter q = (/3 2 — n^k 2 ) 1 / 2 characterizes the decay of 
the field outside the fiber, and the parameter s is de- 
fined as s = (l/q 2 a 2 + l/h 2 a 2 )/[J[(ha)/haJ 1 (ha) + 
K[(qa) I 'qaKx{qa)\, with h = (n 2 k 2 - (3 2 ) 1 ' 2 . The co- 
efficient A is proportional to the amplitude of the field. 
The notation J„ and K n stand for the Bessel functions 
of the first kind and the modified Bessel functions of the 
second kind, respectively. 

In the spherical tensor representation, the components 
of the field envelope are given by 

£-! = V2iA{l + s)K 2 {qry {Pz - 2v \ 
£0 = A^K^qry^-^, 

£ x = -V2iA(l- s)K (qr)e l ' 3z . (14) 

In the case of conventional weakly guiding fibers [2j|, £-\ 
and £q are negligible compared to £\. However, in the 
case of vacuum-clad subwavelength-diameter fibers, £_i 
and £q are not negligible [24[ . 




Radial distance r [|xm] 

FIG. 7: Magnitudes |£-i| (dotted lines), |£o| (dashed lines), 
and (solid lines) of the spherical tensor components of the 
evanescent light field in a fundamental mode with clockwisely 
rotating polarization outside a vacuum-clad subwavelength- 
diameter fiber. The radius of the fiber is a — 0.2 /j,m. The 
light wavelength is A = 685.5 nm (a) and 934.5 nm (b). 

We illustrate in Fig. the magnitudes \£-i\, \£o\, and 
\£i\ of the spherical tensor components of the field out- 
side the fiber. According to the figure, all the three com- 
ponents of the field are comparable to each other in the 
vicinity of the fiber surface. Therefore, we must take into 
account all of these terms when we calculate the light 
shifts of the excited-state levels of atoms outside the thin 
fiber. 

We use eqs. (llUfl and (|14fl together with the Hamilto- 
nian and the hfs energy J2J) to calculate the light shifts 
of the 6P3/2 hfs sublevels of atomic cesium in a two-color 
evanescent field around a vacuum-clad subwavelength- 
diameter fiber. These light shifts are in fact equal to 
the optical potentials of the excited atoms outside the 
fiber. In addition, we also calculate the light shift of the 
6S1/2 state, which is the optical potential of the ground- 
state atoms. To minimize the difference between the light 
shifts of the QP3/2 and 65*1/2 states, we choose the cen- 
tral red- and blue-detuned magic wavelengths for the two 
fields, i.e., Ai = = 934.5 nm and A2 = As = 685.5 
nm. The powers of the laser beams must be chosen ap- 
propriately so that the ground-state optical potential has 
a deep minimum outside the fiber with a high barrier near 
the fiber surface. For this purpose, we choose the pow- 
ers Pi = 11.5 mW and P 2 — 48.5 mW for the red- and 
blue-detuned laser beams, respectively. 

We plot in Fig. E]the light shifts of the 6P 3/2 FM F <-> 
6Si/ 2 F' M' F transitions of the atoms trapped around the 
fiber. The figure shows that the light shifts can be re- 
duced to become less than 14 MHz. Such a shift is com- 
parable to the typical detuning of near-resonant fields 
used for a MOT with cesium atoms (typical detuning is 
10-20 times of the natural linewidth 7 = 5.18 MHz of the 
cesium D 2 line) 0. Thus, due to the use of the red- and 
blue-magic wavelengths, the spatial dependences of the 
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FIG. 8: Spatial dependences of the light shifts of the 
6P 3 / 2 FMf <-> 6S 1/ '2-F , ' M' F transitions of cesium atoms in a 
two-color evanescent field around a vacuum-clad silica-core 
subwavelength-diameter fiber. The radius of the fiber is 
a = 0.2 fim. The two laser beams are tuned to the wave- 
lengths Ai = \r — 934.5 nm and A2 = As = 685.5 nm, 
with the powers Pi = 11.5 mW and P2 = 48.5 mW, respec- 
tively. We show only the results for the transitions involving 
the excited-state sublevels that are split from F = 5. 

light shifts of atomic transitions are weak. This opens up 
an opportunity for state- insensitive trapping, i.e., for si- 
multaneous trapping of ground- and excited-state atoms 
[Hill. If we detune the MOT fields from the cooling 
transition by a negative detuning 5 with a magnitude 
I <5j/2-7r > 14 MHz, then the red-detuning condition for 
the MOT fields (cooling fields) in the presence of the far- 
off-resonance evanescent fields (optical trapping fields) is 
kept throughout the outside of the fiber. This allows the 
simultaneous operation of the MOT and the dipole trap. 

We plot in Fig.[!2]thc optical potentials of the excited- 
and ground-state atoms for the parameters of Fig. [HJ We 
observe from Fig. [5] that the excited- and ground-state 
optical potentials have similar shapes, with deep minima 
located close to each other in space. This indicates the 
possibility of state-insensitive two-color trapping of ce- 
sium atoms around the fiber. Such a state-independent 
trapping scheme allows the simultaneous operation of 
trapping and probing, that is, the operation of trapping 
with continuous observation during the trapping interval. 
Kimble et al. have proposed and demonstrated a similar 
state-insensitive trapping method, which is based on the 
use of counter-propagating laser beams at a red-detuned 
magic wavelength [IJ. 



V. CONCLUSIONS 

In summary, we have shown that the light shifts of 
the 6P3/2-FMF *-* QSi/ 2 F' M' F transitions in atomic ce- 
sium can be minimized by tuning one trapping light 
to around 934.5 nm in wavelength (central red-detuned 
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FIG. 9: Optical potentials (light shifts) of cesium atoms 
in a two-color evanescent field around a vacuum-clad 
subwavelength-diameter fiber. The atoms are in the excited- 
state manifold 6P3/2 ( a ) an d the ground-state manifold 6S1/2 
(b). All the parameters are the same as those for Fig. |H| We 
show only the results for the transitions involving the excited- 
state sublevels that are split from F = 5. 



magic wavelength) and the other light to around 685.5 
nm in wavelength (central blue-detuned magic wave- 
length). We have investigated the light shifts of the ce- 
sium hfs sublevels in a two-color evanescent field around 
a subwavelength-diameter fiber. The simultaneous use 
of the red- and blue-detuned magic wavelengths allows 
state-insensitive two-color trapping and guiding of ce- 
sium atoms along the thin fiber. Our results can be used 
to efficiently load a two-color dipole trap by cesium atoms 
from a magneto-optical trap and to perform continuous 
observations. 
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